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Abstract 

. . . We represent an algorithm reducing a big class of systems of (M + l)-dimensionaI 

i nonlinear partial differential equations (PDEs) to the systems of M-dimensional first 

I order PDEs. Thus, we integrate the original system with respect to only one independent 

CN ■ variable reducing its dimensionality by one. For this reason we call such systems partially 

Ch . integrable ones. In particular, if Af = 1, then M-dimensional PDEs become ODEs. Our 

^ ^ I method may be referred to as a generalization of the method of characteristics. Possible 

■ application in hydrodynamics is discussed. 

00 1 Introduction 

The method of characteristics is known as an effective method for solving the first order scalar 
, multidimensional PDEs reducing them to the appropriate systems of scalar ODEs [Ij. Re- 

garding the systems of scalar multidimensional first order PDEs, their integration is more 
sophisticated. There are few known classes of integrable systems among them [H [21 El IH O [6] . 



> 



Tj- [ We underline only one class of such systems which is most relevant to the algorithm developed 



in our letter, namely, the systems of (M + l)-dimensional PDEs with the same leading parts 



M 



+ 5^a,^ig)-6, =0, j = l,...,iV, (1) 



i=l 



where Oi and bj are arbitrary functions of t, Xj, i = 1, . . . , M and u^^\ J = 1,...,N. Thus, t 
^ ' and Xi, i = 1, . . . , M are independent variables while u^^\ i = 1, . . . , N are scalar fields. It is 
■ shown (see P) that the solutions w*^*^ to the system of nonlinear PDEs ([1]) may be found as 
solutions to the system of non-differential equations 

F('=)(t,x,u) = 0, k = l,...,N, (2) 

where F^'^^ are solutions to the following linear PDE: 

M N 



+ V a-F(^) + V b-F^''^ - 9{F^'^ . . . F^^y) 



Eqs.® 



^ 0. (3) 



1=1 i=l 

We show that the algorithm based on equations ([2]) and ([3]) may be generalized introducing 
a new type of partially integrable systems of multidimensional nonlinear PDEs having the 
following general form: 

M 

Lv + J2 ^^"'^ Lv^r. = A (4) 



m=l 
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where L = dt + YZ=i ^^'"^^x™, v = {v'^^\ is a vector field, i^"^), A and /i^™) (m = 
1, . . . , M) are N x N matrices, A^-dimensional vector and scalars respectively depending on v^^'^ 
and Vx] , i = 1, . . . , N, j = 1, . . . , M. Superscript T means matrix transpose. Here "partially 
integrable systems" means that we are able to reduce the system of [M + l)-dimensional scalar 
nonlinear PDEs to the appropriate system of M-dimensional PDEs, i.e. one can reduce the 
dimensionality of the original nonlinear PDEs by one. In particular, if M = 1, then we reduce 
the two-dimensional system of PDEs to the system of ODEs. Emphasize, that, in general, our 
algorithm provides arbitrary functions of M independent variables in the solution space of 
considered nonlinear PDEs, i.e. we deal with the general solution. This conclusion is justified 
in the examples of Secj3l where the initial value problem for the derived nonlinear systems is 
investigated. 

Among the equations of new type we point on the following (2+l)-dimensional system of 
two PDEs for scalar fields u and w. 

Ut - y^u^i + m^J + ^(m^;, + u^^f + + (m^i - m^Ji + (5) 

+w{i'{u^^ + m^Ji - {u^^ + u^^){u^^ + M^Jtj = 0, 
Iff + 1 = 0, y = const, 

We show that the construction of the general solution to this (2+l)-dimensional system is 
reduced to the solution of the system of two first order two-dimensional nonlinear PDEs: 

H'^^\w + t,u,u^^,u^^,t) + C^^\xi,X2) =0, (6) 

where C'^^^ are arbitrary functions of two variables and the functions H^^"* satisfy the following 
system of linear PDEs: 

Ht\z,uJ'\p^'\t) = (7) 

Hi'\z, uj'\p^^\t) ( - + ^p(^)) + + + (. - tf) , 

H% {z,u, p(^) , p(^) ,t) = Hi'\z,u, p(^) , p(^) ,t)[{u- p(^) - ){z-t) + l), 

H%{z,u, , p^'^ ,t) = H^^\z,u, p(^) , p^'^ , t) - p^^^ - p^'^ ){z-t)-z), 



where k = 1,2. Nonlinear PDEs ([6]) (a) have derivatives with respect to the independent 
variables xi and X2 while variable t appears as a parameter and (6) involve two arbitrary 
functions of variables Xi and X2- As a consequence, the initial value problem may be correctly 
formulated for system ([5]). The peculiarity of system ([5]) is that the properly introduced small 
scale expansion reduces this system to the following single PDE: 

see Sec J3.3[ Eq.® describes, for instance, the special case of the dynamics of two-dimensional 
viscous incompressible fiuid with the constant kinematic viscosity u when two velocity compo- 
nents are equal to each other. Eq.(IH]) becomes (l+l)-dimensional nonlinear Biirgers equation 
(integrable by the Hopf-Cole substitution [3 El [9]) if either Ur^ = or Uy = 0. 

The general algorithm for deriving the new type of systems of nonlinear PDEs is given in 
SecJH Several examples of such systems are represented in Secj3l Concluding remarks are given 
in SecJH 
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2 General algorithm 

2.1 Derivation of partially integrable systems 

Let us consider the (M + l)-dimensional space of variables t and Xi, i = 1, . . . , M, which are 
independent variables of the nonlinear PDEs written for the functions u^'^\ i = 1, . . . , M and 
w^^\ i = 1, . . . , K. Hereafter we will use the notations p'-*-'^ = Ux^, 

u = {m(^) : i = 1, . . . , A^}, w = {w^^^ : / = 1, . . . , K}, x = {xj : j = 1, . . . , M}, (9) 
p={j^^^)^^: z = l,...,iV, ^• = 1,...,M}. 

Theorem 1. Let the functions F^^^ {t, x,u,w,p){k = l,..., M+K) of (M+N+K+NM+l) 
variables satisfy the following overdetermined system of linear PDEs: 

M N K 

Ei.) ,^ Ff) + 5^Fif/.(^-) = J]F(SA« + J]Fl^ii('), (10) 

j=i i=i 1=1 

N K 

j^inmk) F^fi, =^F(gA("™)+^F(^ii(-'), (11) 

i=l 1=1 

which must be compatible. Here h^^\ A^^\ ^(n™)^ ^{nmi) ^ Wi,j,l, n,m, are functions of t, 
X, u, w and p. Let the functions {i = 1, . . . ,N) and w^^'^ (/ = 1, . . . , i^) be solutions to the 
following system of equations: 

F^^\t,x,u,p,w) = Q, k = l,...,N + K, (12) 

where 



^ 0. (13) 

Eqs.(fT2l) 



a(M(i)...«W«;W...wW) 

Then the functions u and w are solutions to the following system of nonlinear (M + 1)- 
dimensional PDEs: 

M N M / M \ 

n=l j=l m=l \ n=l / 

M N M / M \ 

?i=l j = l m=l \ n=l / 

(which is analogous to eq.(jl])). 

Proof. First of all, let us differentiate eq.( fT2il with respect to t and x„, n = 1, . . . , M: 

N N M K 

i=l i=l j=l 1=1 

N N M K 



i=i j=i 
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Consider the following combination of equations ( lT6ll and ( |T71) : 



M 

£W + = 0, k = l,...,N + K. (18) 



the foUowiner linear combination of the derivatives F and F ■ 



Eliminating derivatives and F from system ( fT8l) using eqs. lfTOi) and (fTTi) one obtains 

N / M N M / M \\ 

i=l \ n=l j=l m=l \ n=l / / 

K / M N M / M \ \ 

i=l \ n=l j=l m=l \ n=l / / 

In view of condition f|T3|) . all coefficients ahead of F^a) and F^(i) in eq. f|T9l) must be zero which 
yields the system of nonlinear PDEs fll4lll5p . □ 

Remark: It is obvious, that the system f|T2|) is a system of M-dimensional PDEs, in general, 
due to the functions p'-*-''' = Ux], j = 1, . . . , M. Variable t appears as a parameter in eqs. (fT2|) . 
Thus we reduce the dimensionality of original (M + l)-dimensional nonlinear PDEs fll4fl5p 
by one. It is important that the solution space to nonlinear PDEs does admit the (A^ + K) 
arbitrary functions of M variables in general case. In fact, eqs. fllOflip mean that functions F^''^ 
are arbitrary functions of variables Xi, u^^^ and tw'-'^-* (i = 1, . . . , M, j = 1, . . . , A^, A; = 1, . . . , K). 
In view of condition f|T3|) . one can resolve eqs.f ll2p with respect to u^^^ and w^''^ and rewrite the 
system ([T^ in the following form: 

^0) = F^'^^\t,x,p), j = l,...,N, (20) 
^(fc) _ F(2''=)(t,x,p), k = l,...,K, 

where (A^ + K) functions F^*'-?) keep arbitrary dependence on M variables x^^\ i = 1, . . . , M. 
The availability of arbitrary functions allows one to study, for instance, the initial value problem 
for nonlinear PDEs fll4lll5p . This is reflected in examples of SecEl 

We see that the structure of the system of nonlinear PDEs (1 1411 151) is defined by the coef- 
ficients in eqs. ffTOlITT]) . i.e. by the functions A^'\ A^^"''\ A^^\ i^'""^) and h^^\ These functions 
may not be arbitrary since they must provide the compatibility of the overdetermined system 
ffTOj) and ffTTl) . This problem is discussed in the next subsection. 



2.2 The compatibility conditions of eqs. ( 1101) and ( 1111) 



The compatibility conditions of eqs. (|T0|) and (|TT1) read: i^^^"™-^) = F^^]„), n = 1,...,N, 



m = l,...,M, k = l,...,N + K,oT 

N K M 



E + E ^^"'"'^ + E F^^'^ = 0, (21) 

i=l 1=1 j=l 
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where 



Q{nmi) ^ ^ 



(i) 



,(nm) 



(nmi) 



M N 

- E + E {^''^' - A$,A^) + (22) 

i=i i=i 



nml) 



1=1 



M 



N 



_ A 

,(nm) 



{nml) 



i=\ 



Q{nml) ^ ^(0 
[=1 

AT K 



E 47'^ /i^^'^ + E - ^S)^^""'^) + (23) 

i=i 



(24) 



If at least some coefficients are nonzero in eq.f l2T|) . then this equation exhibits a hnear depen- 
dence among the derivatives F^{i), F^{i) and F^ , V i, j, /. Since derivatives F^(i) and F^n) are 
hnearly independent (in view of condition f lT^ ). the above hnear dependence would reduce the 
dimensionality of nonlinear PDEs, which is not desirable. Thus it is reasonable to require that 
all coefficients in eq. fl2T]) are zeros: 



Q{nmi) Q{n'ml) g 

g(«™j) = 0, z,n = l,...,N, j,m = l,...,M, 1 = 1,. ..,K. 



(25) 
(26) 



Equations fl25|26p represent the ffist system of nonlinear PDEs defining the functions A^'^\ A^^\^ 



Next, the compatibility condition of eqs.f lTT]) reads: E^^ 
s,m = 1, ... ,M, k = 1,. .. ,N + K, or 



(rsk) 



nm) 



j-,(nmk) -I 

K(rs) % r,n = 1, . . . ,iV, 



N 



K 



Ep{k) g(rsnmi) , \ ^ jpi^) ^(rsnml) _ g 



(27) 



where 



(=1 



N 



p{nm) p{rs) / j ^ ^{i) ^(i) J 

i=l 

K 

[=1 

N 

oirsnml) _ A{rsl) A{nml) ST^ f A(nml) .(rsi) A^^^^-) A(nmi)\ ^ 



(28) 



(29) 



i=l 



K 



E( J^inml) j^{rsl) _ ^M)^C 



nml) 



1=1 
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Similar to eq. (!2T]) . eq. (!27|) represents the linear relation among derivatives F^^]^ and F^^j^ (V 
i,l,k). Again, condition (fT3|) suggests us to take all coefficients in eq.( !27|) equal to zero: 

g{rsnmi) ^{rsnml) g (30) 

i,r,n = 1, . . . , N, s,m = 1, . . . , M, l = l,...,K. 

Equations (130|) represent the second system of nonlinear PDEs defining the functions and 
^{nmi)^ Thus, if the functions A^'\ A^^\ j^inmi)^ ^{nmi) g^^^ ^{j) satisfy the system of nonlinear 
PDEs fl25|26ll30l) . then homogeneous system of linear PDEs (fTOl) and (fTTl) is compatible. 

2.2.1 Integration of system (1251126030!) 

Hereafter we consider such potentials in eqs.f llO|[TT!) which do not depend on x explicitly. Then 
the terms with Xj-derivatives disappear from eqs. fl22ll23p . In order to write the system of 
nonlinear PDEs fll4lll5p in the explicit form we have to find the admissible expressions for the 
coefficients in the linear system fllOlllip . i.e. we have to solve the system fl25ll26|30p . which is 
nonlinear in general. In the simplest cases (see Sec J3.3p . this system reduces to the linear one 
with constant coefficients which may be integrated directly. 

In general, the following theorem provides us a big manifold of acceptable functions A^'^\ 
A^^\ A^'^™) and 

Theorem 2. Let \E'(^\ = 1, . . . , N + M, be arbitrary functions of variables t, m, w such 
that the system 

N K 

ffik) _ A^^''^ = 5^I^„w*^''^^^'^ + 5^^^(0^^''^^^'\ (31) 

i=l 1=1 

N K 

i=l 1=1 

n = l,...,N, m = l,...,M, k = l,...,N + K 
can be viewed as a uniquely solvable system of algebraic equations for the functions A^^\ A^''\ 



^ 0. (33) 



Here we introduce operators 

= -i^ik) ^ ^(y;k)^{k) ^ (34) 

where y is one of the variables t, w'^''\ p^'^^\ Then functions A^*^ and 
satisfy the system (I25|30p . 

Proof. To prove this theorem one has to show that system fl25|30p is the compatibility 
condition of the overdetermined system f l3T][32|) . In fact, the compatibility condition of eqs.f l3ip 
and ([32]) reads: Dp(™) = V n, m, fc, or 

N K 

E ^SQo"™^ + E = 0' (35) 

i=l 1=1 

n = l,...,N, m = l,...,M, k = 1,...,N + K, 
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where Qq"™-* and Qq"™''' equal to and ^(""^0 respectively without terms having deriva- 

tives with respect to x, see eqs. fl22|23p . 

Next, the compatibility conditions of eqs. (Eg) read: = V n, m, r, s, /c, 

or 

JV K 

J2 ^Iw 5^™^) + J2 ^S) ^^"'"""'^ = 0, (36) 

i=l 1=1 

r,n = l,...,N, s,m = l,...,M, k = l,...,N + K, 

where coefficients 5'(^«"™) and 5'(^*"-'"0 are give by eqs.( l28|29l) . Both eqs.fl35l) and ( l36l) may be 
viewed as linear homogeneous systems for the functions Q^q^^\ Qq"™'"* and S^^^"'"^^^ 
respectively. However, due to condition ( l33ll . the homogeneous systems (l35ll and ( l36l) have only 
zero solutions, which yield the systems (I25|30p . □ 

In turn, the functions h^^\ j = 1, . . . , M, may be found as solutions to the linear PDEs fl26|l . 

Remark: Operators Dy are introduced in order to provide us the possibility to remove a 
dependency on some particular variable. For instance, if we would like to eliminate explicit 
dependency of functions A^'\ A^^\ ^(nmi)^ ^{nmi) ^.^^ /^(j) 

on some particular variable yo, then 
we must remove derivative in the appropriate operator Dyg which acquires the following form: 

2.2.2 Admissible reductions for the potentials of eqs. (1101) and (111 ft 

The described method for construction of potentials in eqs. (fTO|) and (fTT]) preserves (A^ + K) 
arbitrary functions '^^'^^ of all variables t, u, w, p. The availability of such freedom suggests 
us to look for additional compatible relations among the above potentials. A big class of these 
relations may be established by means of the spectral theory developed for the nonlinear PDEs 
assotiated with commuting vector fields |TOl|TT]. Namely, the following theorem is valid. 
Theorem 3. Let us consider the spectral problem for some spectral function x(A): 

N K 

T^W := x*(A) + Ax,oo(A) = $^I^„«x(A)A« + 5^D^(ox(A)i('), (37) 

i=l 1=1 

N 

^^"'"^ := Xp(""^)(A) + Ax,_(A) = 5^D,„x(A)A("™)+ (38) 

i=l 

K 

5^D^(0X(A)i(""^^ n = l,...,A, m = l,...,M, 

1=1 

where rjij are additional parameters, which, in particular, may coincide with some parameters 
t, M, w or p; A is a complex spectral parameter. The spectral function x depends on A as well as 
on all other parameters such as t, it, p, w and r]ij. Then functions A^^ , A^^\ ^{nmi) ^ ^(nmi) 
h^^^ satisfy the system ( 125|) (without terms having Xj-derivatives) and (!30|) as well as additional 
reductions 

j^ii) _ j^inmi) ^(l) _ ^{nml) j^{rsi) _ ^(nmi) ^{rsl) _ ^(nml) /ggN 

n,r,i = 1, . . . , N, s, m = 1, . . . , M, I = 1, . . . , K. 
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Proof. Overdetermined linear system ( l371l38l) for the function x must be compatible. Com- 



patibility conditions of eqs.(|371l and ( l38ll read: 1-6%^^.) = 'H^t'^\ V n, m, or 



N K 

E D^^.X^'^Qt"^'^ + ^^(oX^'^Qr'^ + (40) 

i=l 1=1 

(N K \ 

i=l 1=1 J 

n = l,...,N, m = l,...,M, k = 1,...,N + K. 
Next, the compatibility conditions of eqs.( 15K]l read: Dp(nm,)H^'^^'^ = Dp(rs)H^""^\ yn,m,r,s, or 

N K 

E + E ^«;(0X^''^'5^"""'^ + (41) 

1=1 1=1 

/ N K \ 

\i=l 1=1 ) 

r,n = l,...,N, s,m=l,...,M, k = l,...,N + K. 

Following the general strategy of spectral theory for the equations associated with commuting 
vector fields [TOl |TT] we conclude that all coefficients in eqs. fl40|4ip are zeros, which yield 
relations along with eqs. fl25]l30|) . □ 

Again, the appropriate functions h^^'^ may be found as solutions to the linear system f l2B]) . 

Remark: Relations f l5I?]) may be satisfied, if we introduce functions A^^^ and A^^'^ by the 
following formulas: 



3 Simple examples of nonlinear PDEs and implicit rep- 
resentation of their solutions 

After functions A^^ , A^^\ ^('^™), i('^™0 and /i^^^ have been constructed, one has to find {N + K) 
solutions F'^^\ k = 1, . . . ,N + K, to the linear system ( ITOlITT]) such that condition (fT3ll is 
satisfied. This problem may be solved numerically, in general. The last step in construction 
of solutions to nonlinear PDEs fll4lll5p is solving the system of M-dimensional PDEs ( IT2l) for 
the functions u and w. Some examples of nonlinear PDEs together with implicit description of 
their solutions are represented below. 



3.1 Example 1: the system of first order nonlinear PDEs with the 
same leading part 

We show that our algorithm reduces to the method of characteristics in the particular case. 
Let = = and F^''^ be independent on p. Then eq.l ITT]) disappears, so that the 

functions F^^^ become solutions to the single linear PDE (fTOj) : 

M N K 

Ft' + E ^if^'' = E ^S^'^ + E )^'^ k = l,...,N + K (43) 
j=i i=i 1=1 



with arbitrary functions A^'^\ A^^^ and h^^^ of all variables t, x, u and w. Eq.(l43l) is equivalent 
to eq.(|3]). In turn, nonlinear PDEs ( !T4fT5|) read 



M 



+ + =0, (44) 

n=l 
M 

'^'^ + $^«^2^^"^+^^'^ =0. (45) 



n=l 



This is a system of nonlinear PDEs with the same leading part, which is analogous to system 
([1]). We do not discuss its solutions. 

3.2 Example 2: nonlinear PDEs in the form + f^'^\u,p)^ = 0, 

In this subsection we consider an example of the system of nonlinear PDEs which may be given 
the following general form: 

t + (m,p))^ = 0, k = l,...,N. (46) 

Of course, this system may be immediately integrated over t yielding t + f^''\u,p) = c^'^^x), 
where c^^-* are arbitrary functions of x, so that we do not need to apply our algorithm. However, 
we represent derivation of such equations using our algorithm in order to demonstrate its 
relation with the classical integration methods. 

Let K = and be the following linear functions of u and p (independent on w): 

N M 

^(k) ^ ^ik) + Y,Y1 A; = 1, . . . , A^, , a^'^'^) = const. (47) 

i=i j=i 

We take 

= A^'^^vi/W, A("'"') = 0, A(^^"™^^'^) = A("™'=), /i^™) = 0, (48) 
k,n = 1, . . . ,N, m = 1, . . . ,M. 

Algebraic system fl31|32p becomes simpler in this case: 

N 

X(>')^W = J2^[%A^^, (49) 

1=1 

N 

D,,_,v1>W = ^v1>Wa("™), k,n = l,...,N, m = l,...,M. (50) 

i=l 

Its solution reads: 

^(fc) = A(fc)^(fc)^ Ainmk) ^ ^{nmk) ^ ^{nmk)^{k) ^ k,n=l,...,N, m = 1, . . . , M. (51) 

Regarding the nonlinear system, we observe that eq. ffTSj) disappears while eq. ffT^ reads: 

N M 

^« + + ^ ^ 0, k,n = l,...,N, m = l,...,M. (52) 

j=l m=l 



Now we consider the system ( |TO|TT]) defining functions F^''\ k = 1, . . . , N . Introducing new 
variables instead of u^^\ i = 1, . . . , N, we rewrite this system as follows: 

N N 

f}'^ = E ^-S^^^^^^ ^il) = E i^£^^^A(— ), k,n = l,...,N, m = l,...,M. (53) 

i=l i=l 

Solution to system (!53|) reads: 

N M 

F(*^) = G(^)(xi,...,a;M,ln^(^) + A«f + EE^'^''^^^'''^'---' ^^^) 

i=i j=i 

N M 

In v^(^) + A(^)t + E E I'^'^'^A^^^-^)) , k = l,...,N, 
i=i j=i 

where G*^^-' are arbitrary functions of M + variables. We may choose functions G^'^^ such that 
the system G^'^^xi, . . . , xm, -^i, • • • , z^) = 0, k = 1, . . . , N, (see eqs. (fT^ ) is uniquely solvable 
with respect to variables Zi, . . . , z^. Then system fll2p may be replaced with the following one: 

N M 

ln^W + A('=)t + EE^^'^^^^'^''^ = '^^''^(^i'---'^^^)' k = l,...,N, (55) 
i=i j=i 

where C^*^ are arbitrary functions of M variables. It is obvious that the t-derivatives of eqs.f l55p 
yield the system ( 152|) . so that we deal with the nonlinear PDEs writable in form (H6|l . Eqs.f l55p 
are M-dimensional nonlinear PDEs where t appears as a parameter. 

The above arbitrary functions C*-*^-* may be used in order to solve the initial value problem. 
In fact let Uq = u\t=o be the initial condition. The appropriate initial values p^-^^ and \l/o^ for 
the functions p^*-'-' = Ux] and may be calculated directly. Putting t = in eqs. (l55|) one gets 
the following equations: 

N M 

In^i'^) + ^^pfe)Afe-^) = C'^'\x„ . . . ,xm), k = l,...,N, (56) 

i=i j=i 

which uniquely define functions C^^\ k = 1, . . . , N. 

In particular, let N = 1, M = 2, u = m^^), A = A^^), so that ^^^^ = u + a'^'^^'^^^i + a^^^^^w^a 
and eq.( l52|) reduces to the single nonlinear PDE: 

ut + X{u^^ +u) + u^^t + {u^i + u)u^2t = 0, (57) 

where we introduce the new variables (^j, i = 1,2 in accordance with the following formulas: 
= + a^^^^^dx^, d^^ = X^^^^^d^, + X^^^^^d^,. The system §Si> reduces to the single 

two-dimensional PDE, where t appears as a parameter: 

\n{u^, +u) + Xt + u^, = C(i)(6,6)- (58) 

If there is no dependency on ^i, i.e. u^-^ = 0, then eq. flFH) reads 

Ut + uu^^t + Xu = 0. (59) 

In this case eq.([58D becomes ODE, C(^)(^i,6) ^ <^(6): 

lnu + Xt + u^^ = Ci^2), (60) 

where C(^2) is an arbitrary function of single variable. 
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3.3 Example 3: a new class of nonlinear PDEs 

In this section we consider an example when equations (125!) . (!26|) and (130!) reduce to the hnear 
ones with constant coefficients and consequently may be easily integrated. In addition, we 
assume that all functions appearing in system of equations ( 125|) . fl26!l and (130|) do not depend 
on t, Xj and Ui, so that the appropriate derivatives disappear from these equations. 

3.3.1 The general solution to eqs. and (ED]) with i^*) = const and i^"*"^) = 

const and appropriate nonlinear PDEs. 

Let = 7^"™) = const and i^^) = 7^^^ = const. Then equations ([25]), ([26]) and ([30]) reduce 

to the following system: 

^^ju + x: (a^-^s^') - = 0, (61) 

^p(U-E^I:i)7("'"'^ = 0, (62) 

A^r2) - + E «r^7^^^') - A^i^i"^') = 0. (63) 
[=1 

Solution to eq.( l62|) reads: 

h^'^=g^'HO, j = l,...,K, (64) 

{N M 
^(i) = ^{i) + T^"™^^""'^ ■.i = l,...,K 

n=l m=l 

System fl61f63p may be readily integrated. Assuming that (at least) some 7*^*^ are nonzero one 
obtains: 

/(K N M \ 

^i=^«-^ ^+^^(0, (65) 

N M \ 

E E E k<.nh)p(Hn) 

J '1=^ '1=^^!=^ J +G^^'"'\ (66) 

a«(A;) I ^ A;«7("'"') - k^^'^n 
a("'"')(A;) = ^^^^i^ (67) 

1=1 

i,n = 1, . . . , N, m = 1, . . . , M, 

A; = {A;(^) : z = 1, . . . , A:^"'") : n = 1, . . . , iV, m = 1, . . . , M}, 



where a; (A;) is some measure and F is a domain in the space of vector parameter k; G\ 
const., G^'^\^) are arbitrary functions of K variables ^'^'^\ i = 1, . . . , K. To satisfy eq. fl67l) we 



{nmi) 
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take the following structure of functions a^"''^^\k) and a^^\k): 



K 



, 1=1 



K 



'\k) = h^\k)Y^k^'^-i^'\ (69) 
1=1 

where b^^\k), i = 1, . . . , N, are arbitrary functions. 

In particular, if all j^^^ = 0, then a^''\k) = and eq. fl65l) reduces to the following one: 

A« = ^^(e), (70) 

while eq. fl67|) must be replaced with the following relations among a("-"^*): 

which agrees with eqs.( l68l) . Nonlinear system ( |T4l[T5l) reads as follows: 

M N M / M \ 

^ + E ^9^(0 + ^« + E E ^^^'^^ + E = 0, (72) 

n=l j=l 771=1 \ n=l / 

M 

+ E^2^?^"^(^)+^^'^ = 0' ^ = i.---,^^: / = i,...,ir. (73) 

n=l 

Finally, eqs. (ITnllTT]) read: 

M AT K 

Ft' + E ^if ^"'^ = E ^iS^'^ + E ^(i^'^ (74) 



j=i i=i /=i 

AT K 



F^L = E ^SU^"™^ + E fI%7^-'\ (75) 

j=i 1=1 
n = l,...,N, m = l,...,M, k = 1,...,N + K. 



Solutions F^''^ to system fl74ll75p must be used in eqs. (fT2!) for construction the functions u and 
w. 

3.3.2 An application to the hydrodynamics of viscous fluid 



Now we consider a simple particular solution to the system fl61tl63p together with the appropriate 
nonlinear PDEs in more details. This case agrees with the general solution derived in Sec J3. 3. 11 
Let us take N = K = 1, M = 2, u^^^ = u, w^^^ = w, p^"*") = = u^^, h^^ = 0. Then 
eqs. (l6m63|) reduce to the following ones: 

+ AL^'-^W - A^^V'"-'' =0, m = 1, 2, (76) 
- <\f^ + A^^-^'h^'-'^ - AL^^S^^"^) =0, m = 1, 2. (77) 
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In turn, the system of nonlinear equations ( 172f73l) gets the following form: 

2 

m,=l 
2 

y^t + 7^'^ + E 7^'™'^^.™* = 0. (79) 

m=l 

The general solution to eqs. fl76|77p reads: 

A('-) = --La;M«, + 2-^4('> + cW, (80) 

where A^^-* is an arbitrary function of t/;, p*^^-* and p*^^^ while c^*"^ (m = 1,2) are arbitrary 
constants. We look for the solution to linear system (I74II75P in the following form: 

2 

p(k) ^ i/W(^^^^p(i)^p(2)^^) + C'('=)(a^i,a;2), z = 1(7 + ^p('")7(^'"^) +7«t, (81) 

m=l 

where C'^^^ (A; = 1,2) are arbitrary functions of two variables while H'^^^ (/c = 1,2) satisfy the 
following linear system: 

hP = ffW (a(i)L=,_j^^^^p(„)^(.™.)_^(i),) , (82) 
H% = H^'^ , m = 1,2. (83) 

It is simple to show that, if A^^^ does not depend on w, then system fl82|83p may be integrated 
analytically and the t-derivatives may be removed from nonlinear PDEs ( 178)1791) by the simple 
method discussed in Sec J3.2[ In general, one has to solve system ( 182|83|) numerically, so that 
the appropriate nonlinear PDEs are beyond the family of PDEs considered in Sec J3.2[ Namely 
this case leads to the new partially integrable nonlinear PDEs. 

The final step is solving the two-dimensional first order PDE (|T2l) . which reads as follows: 



H'^''\z,u,p'^^\p'^^\t) + C'^''\xi,X2) = 0, A; = 1,2. (84) 
In addition, functions F^''^ must satisfy the condition f|T3l) which now reads: 

9(i7«,i/(2)) 



^ 0. (85) 

Eqs.( 



d{z, u) 

Without loss of generality, we take H^'^^z, u,p^^\p^'^\ t) = z, so that system (|8^ can be written 
as follows: 

H^'\z,u,p^'\p^^\t) + C^'\x^,X2) = 0, (86) 

2 

z + C(2)(xi, X2) = w + Y, p^'^h^^""^^ + l^^'^t + C(2)(xi, X2) = (87) 



m=l 



It is obvious that the described algorithm preserves two arbitrary functions C^'^\xi,X2) (m = 
1,2) of two variables in the solution space to nonlinear system ( |78|79l) . These functions may 
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be uniquely fixed by the initial conditions at t = as follows {uq = u\t=o, wq = w\t=o, Po^^ = 

2 

m=l 

Thus, our algorithm is suitable for the initial value problem. 

Now we turn to the more specific example. Let us take A^^'^ in the following form: 

= + + 1(^(1) + p(2))2 ^ ^3 (89) 

and use the following values for the constant parameters: 

c(i) = l, c(2) = -z, 7« = 1, ^(111) = ^(121) = 0. (90) 
Then eq. (l80|) yields the following expressions for A^^'^^\ m = 1,2: 

A(lll) = {U- - + C(l), yl(121) ^ _ _ ^(2))^ ^ J2)_ (g^) 

Substituting eqs.( l89ll9Ti) into eqs.l lTHi) and ( 1791) one obtains the following nonlinear system: 

Ut - 1^{Uxr + ^Wxa) + ^(Mxi + M^a)^ + + {u^^ - iUx^)t + (92) 



0, 



+ 1 = 0. (93) 



In order to construct the general solution to system fl92|93l) . one has to substitute eqs.( l89ll9T]l 
into the system fl82|83l) (see eqs.©) and solve it for H'^^\ = 1, 2. After that, the initial value 
problem for eqs. (]92|93p may be solved by means of eqs.( j86ll88l) . System (1921193^ is represented 
in the Introduction, see eqs. ([5]). 



Small scale expansion of eqs. (l921l93l) and (2+l)-dimensional dynamics of viscous 

fluid. Let us apply the differential operator [dx^ + dx2) to the nonlinear PDE ( l92l) and intro- 
duce a small parameter e by the following formulas: 

Ux^ + Ux2 = £v, Vx, = ev^^, w = e^w, Vt = evr,, (94) 
where rj, C,i and ^2 are the new independent variables and v, w are the new fields. One obtains 
e^(v^ - u{v^^ + w^jj + e^(v{v^, + v^J + v^.r, - iv^2ri) = o{e^). (95) 

We see that the field w does not appear in two leading terms of eq. fl95|) . 

Now we eliminate r^-derivative in the last term of the RHS of eq. fl95l) substituting v^j = 
y{v^^ + zf^J + •^(^)- Finally, let us introduce the independent variables r, yi and y2 instead of 
77, (^1 and ^2 by means of the following relations: — i^(f^i + iv^^^) = evr, v^^ = Vy-, i = 1,2. 
Then the leading term of eq.f l95l) yields the following nonlinear PDE for the function v: 

Vr + V{Vy, + Vy^) - iy{Vy^y, + V y^y^) = 0, (96) 
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which describes, for instance, the two-dimensional dynamics of viscous fluid with equal com- 
ponents of velocity. Here i/ is the constant kinematic viscosity. This equation is written in the 
Introduction, see eq.®. 

Remark 1: The small parameter must be introduced not only in eq. fl92!) . but also in 
eq.( p3|) . Introducing e in the solution to eq. (l93|) (see eq. (!87|) with constant parameters given 
by formulas (MI]) ) we obtain 



w 



CnUyi_^^y_l_^^xr (97) 



This formula requires the special form of the arbitrary function C-^-*. But, as far as the field 
w does not appear in eq. fl95|) as well as in the final nonlinear PDE fl96|) . we do not consider 
eq.f l97|) in more details. 

Remark 2: It is important to note that equation fl96|) may be derived as a small scale 
expansion of a big family of more complicated nonlinear PDEs inside of the class of partially 
integrable PDEs. The system fl92ll93p is a representative of this family produced by our choice 
of the function AS^^ and constant parameters, see eqs. (l89|l and (!90l) . Another representative is 
produced by the following choice of the function A!>^^: 



w\ / 



with constant parameters given by formulas (1901) . Thus, eq.( l96l) is, in some sense, "an attractor" 
in a big family of partially integrable nonlinear PDEs. 



4 Conclusions 

We suggest an algorithm allowing one to reduce a big class of systems of {M + l)-dimensional 
scalar nonlinear PDEs to the appropriate systems of nonlinear M-dimensional PDEs. In partic- 
ular, if M = 1, then the M-dimensional system becomes the system of ODEs. We refer to these 
PDEs as partially integrable ones since we are not able to integrate them over all variables. 
Among examples, there is a family of (2+l)-dimensional PDEs having the same small scale 
hmit applicable to the hydrodynamics, see eq.f l96|) . Existence of this example suggests us to 
look for further improvements of the represented algorithm in order to integrate eq. (l96p directly 
and exhibit new nonlinear PDEs which are both applicable and integrable. The represented 
method may be viewed as a generalization of the method of characteristics. 

This work is supported by the RFBR grants 10-01-00787 and 09-01-92439 and by the grant 
NS-6885.2010.2. 
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